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Bacteroides timonensis strain API^ (= CSUR PI 94 = DSM 26083) is the type strain of B. 
timonensis sp. nov. This strain, whose genome is described here, was isolated from the 
fecal flora of a 21 -year-old French Caucasoid female who suffered from severe anorexia 
nervosa. Bacteroides timonensis is a Gram-negative, obligate anaerobic bacillus. FHere we 
describe the features of this organism, together with the complete genome sequence and 
annotation. The 7,130,768 bp long genome (1 chromosome, no plasmid) exhibits a Gh-C 
content of 43.3% and contains 5,786 protein-coding and 59 RNA genes, including 2 
rRNA genes. 



Introduction 



Bacteroides timonensis strain APIt [= CSUR P194 
= DSM 26083) is the type strain of B. timonensis 
sp. nov. This bacterium was isolated from the 
stool sample of a 21-year-old French Caucasoid 
female in an effort of cultivating individually all 
bacterial species within human feces [1]. It is a 
Gram-negative, anaerobic, indole-positive rod- 
shaped bacillus. 



support the circumscription of the type species, 
B. timonensis. 



The conventional genetic parameters used in the 
delineation of bacterial species include 16S 
rRNA sequence identity and phylogeny [2,3], ge- 
nomic G + C content diversity and DNA-DNA hy- 
bridization [DDH) [4,5]. These tools have Hmita- 
tions, notably because their cutoff values vary 
across species or genera [6]. With the introduc- 
tion of high-throughput sequencing techniques 
[7], a wealth of genomic data was made available 
for many bacterial species. We recently pro- 
posed to include genomic data in a polyphasic 
approach to describe new bacterial taxa 
[taxono-genomics) [8]. This strategy combines 
phenotypic characteristics, notably the MALDI- 
TOF MS spectrum, and genomic analysis [8-37]. 

Here, we present a summary classification and a 
set of features for B. timonensis sp. nov. strain 
APIT (;= CSUR P194 = DSM 26083) together with 
the description of the complete genome se- 
quencing and annotation. These characteristics 



The genus Bacteroides (Castellani and Chalmers 
1919) was created in 1919 [38]. Currently, it is 
one of the largest genera among the human gut 
microbiota [39], and consists of 91 species and 5 
subspecies with vaUdly pubUshed names [40]. 
Bacteroides species are Gram-negative, non- 
spore-forming, non-motile and anaerobic rods 
that are generally isolated from the gastrointes- 
tinal tract of mammals [41]. They have symbi- 
otic relationships with humans and play many 
beneficial roles on normal intestinal physiology 
and function. Several Bacteroides species are 
identified as opportunistic pathogens when iso- 
lated from anaerobic infections [42]. 



Classification and features 



A stool sample was collected from 21-year-old 
French Caucasoid female who suffered from se- 
vere restrictive anorexia nervosa from the age of 
12 years. At the time of sample collection, she 
had been hospitalized for recent aggravation of 
her medical condition [BMI: 10.4 kg/m^). The 
patient's written consent and the agreement of 
the local ethics committee of the IFR48 [Mar- 
seille, France) were obtained under agreement 
number 09-022. The feces sample of this patient 
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was stored at -80°C immediately after collection. 
Strain APIt (Table 1) was isolated in November 
2011 after 1 month of incubation in Columbia 
agar (BioMerieux, Marcy I'Etoile, France]. Sever- 
al other new bacterial species were isolated 
fr-om this stool specimen using various culture 
conditions. 

When compared to sequences available in 
GenBank, the 16S rRNA gene sequence of B. 
timonensis strain APli" (GenBank accession 
number JX041639) exhibited an identity of 



97.00% with Bacteroides cellulosilyticus (Figure 
1). This value was the highest similarity ob- 
served, but was lower than the 97.8% 16S rRNA 
gene sequence threshold recommended by 
Stackebrandt and Ebers (2006) to delineate a 
new species without carrying out DNA-DNA hy- 
bridization [3], and was in the 74. 8 to 98.7% 
range of 16S rRNA identity values observed 
among 41 Bacteroides species with validly pub- 
Hshed names [56]. 



Table 1. Classification and general features of Bacteroides timonensis strain API^ according to the MIGS 



recommen 


dations [43] 






MIGS ID 


Property 


Term 


Evidence code' 




Current classification 


Domain Bacteria 


TAS [44] 






Phylum Bacteroidetes 


IAS [45,46] 






Class Bacteroidia 


TAS [45,47] 






Order Bacternidale^ 


TAS [45,48] 






Familv Ri^rfproiHi^cPi^p 

1 ci 1 1 1 1 1 y c-4 cv^ 1 v_/ / cx 


TAS 149 501 






("iPnim R^rtprnirJp^ 


IDA 149 51-541 






Species Bacteroides timonensis 


IPl A 






Type strain API^ 


IDA 




Gram stain 


Negative 


IDA 




Cell shape 


Rod 


IDA 




Moti 1 itv 


Non mntilp 

L 1 \-/ 1 1 III \-/ L 1 1 V . 


IDA 




Sporulation 


Non sporulating 


lUA 




Temperature range 


Mesophile 


IDA 




Optimum tempera- 








ture 


37°C 


IDA 


MIGS-6.3 


Salinity 


Unknown 


IDA 


MIGS-22 


Oxygen requirement 


Anaerobic 


IDA 




Carbon source 


Unknown 


IDA 




Energy source 


Unknown 


IDA 


MlCS-6 


Habitat 


Human gut 


IDA 


MlCS-15 


Biotic relationship 


Free living 


IDA 




Pathogenicity 


Unknown 






Biosafety level 


2 




MlCS-14 


Isolation 


Human feces 




MlCS-4 


Geographic location 


France 


IDA 




Sample collection 






MlGS-5 


time 


November 201 1 


IDA 


MlGS-4.1 


Latitude 


43.296482 


IDA 


MlGS-4.1 


Longitude 


5.36978 


IDA 


MIGS-4.3 


Depth 


surface 


IDA 


MIGS-4.4 


Altitude 


0 m above sea level 


IDA 



Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report 
exists in the literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, iso- 
lated sample, but based on a generally accepted property for the species, or anecdotal evidence). These 
evidence codes are from the Gene Ontology project [55]. If the evidence is IDA, then the property was di- 
rectly observed for a live isolate by one of the authors or an expert mentioned in the acknowledgements. 



1182 



Standards in Genomic Sciences 



Ramasamy et al. 




Bacteroides ceUulosilyticus (A| SB 32 43] 
Bacteroides intt^stinalis [AB2 143ZB] 
Bocternides oleiciplonua 

Bacteroides helcvgenes (CP002352] 

Bacteroides Tiinofterisis (JX041 639) 

Bacteroides /luxus {AE'i'iGSOl) 

Bacteroides unifcrtnis (AE0501 10] 
BacKmides eqijerthii (AROSfi 107) 
Bacteroides clams [AB490801) 
Bacteroides stercoris 
Bacteroides aiccae (XE39S1] 

Bacteniidts wlyersiue (A¥60eS96) 
Ractuwides nnrdsi ( AY608697) 
Bacteroides fragiUs [CR626927) 



Bacteroides ovatus (ABOSOlOfiJ 
Bacteroides xylanisolvens (AM Z3 0650] 

Bacteroklei- Jin egaldii [AB2 Z2 699) 
Bacteroides thetaictaomicron [AEO 15923) 
Bacteroides faecis [GQ496624) 

Bacteroides coprosuis (AF3i9778) 

Bacteroides vulgatus (CP000139) 

Bacteroides dorei {AB242 14-2) 
Bacteroides wassiliensis [ABS 1 Q703) 
Bacteroides roproco to [ A B 2 0 02 2 4) 
Bacteroides plelieiiis (AB2t)0217) 
Bacteroides saldttitmnis [CP002&30) 



Bacteniides ctiprophiSus (AB2S0026) 

Prevotello melaninogenica (AY323525) 



Figure 1. Phylogenetic tree highlighting the position of Bacteroides timonensis strain APl^ rela- 
tive to other type strains within the Bacteroides genus. GenBank accession numbers are indicat- 
ed in parentheses. Sequences were aligned using CLUSTALW, and phylogenetic inferences 
were obtained using the maximum-likelihood method within the MEGA software. Numbers at 
the nodes are percentages of bootstrap values obtained from 500 replicates. Prevotella 
melaninogenica was used as outgroup. The scale bar represents a 2% nucleotide sequence di- 
vergence. 



Four different growth temperatures [25, 30, 37, 
45°C) were tested; growth occurred between 25 
and 37°C, but optimal growth was observed at 
37°C, 24 hours after inoculation. No growth oc- 
curred at 45°C. Colonies were translucent and 
approximately 0.3 mm in diameter on 5% sheep 
blood-enriched Columbia agar [BioMerieux). 
Growth of the strain was tested in the same agar 
under anaerobic and microaerophihc conditions 
using GENbag anaer and GENbag microaer sys- 
tems, respectively [BioMerieux), and under aer- 
obic conditions, with or without 5% CO2. Growth 
was observed under anaerobic and 
microaerophilic conditions, and only weakly 
with 5% CO2. No growth occurred under aerobic 
condition without CO2. Gram staining showed 
short Gram-negative rods unable to form spores 
[Figure 1). A motility test was negative. Cells 
grown on agar are translucent and exhibit a 
mean diameter of 0.88 |im in electron microsco- 
py [Figure 2, Figure 3). 

Strain APIt exhibited catalase but no oxidase ac- 
tivity [Table 2). Using an API Rapid ID 32A strip 
[BioMerieux), positive reactions were obtained 



for arginine dihydrolase, a-galactosidase, (B- 
galactosidase, a-glucosidase, (B-glucosidase, a- 
arabinosidase, N-acetyl-(B-glucosaminidase, glu- 
tamic acid decarboxylase, a-fucosidase, nitrate 
reduction, indole production, alkaline phospha- 
tase, proline arylamidase, leucyl glycine 
arylamidase, alanine arylamidase, glutamyl glu- 
tamic acid arylamidase, and fermentation of 
mannose and raffmose. Weak activities were ob- 
served for glycine arylamidase and serine 
arylamidase. Negative reactions were obtained 
for urease, (B-galctosidase-6-phosphatase, (B- 
glucuronidase, arginine arylamidase, phenylala- 
nine arylamidase, leucine arylamidase, 
pyroglutamic acid arylamidase, tyrosine 
arylamidase and histidine arylamidase. Using an 
API 50CH strip [Biomerieux), strain APIt was 
asaccharolytic. B. timonensis is susceptible to 
amoxicillin-clavulanate, ceftriaxon, imipenem, 
trimethoprim-sulfamethoxazole, metronidazole 
and doxycycline but resistant to amoxicillin, 
vancomycin and gentamicin. By comparison 
with other Bacteroides species, B. timonensis dif- 
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fered in production of indole, nitrate reductase, (B-galactosidase and acidification of sugars. 




Figure 2. Gram staining of B. timonensis strain API^ 




Figure 3. Transmission electron microscopy of B. timonensis 
strain AP1\ made using a Morgani 268D (Philips) at an operat- 
ing voltage of 60kV.The scale bar represents 200 nm 
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Matrix-assisted laser-desorption/ionization 
time-of-flight (MALDI-TOF) MS protein analysis 
was carried out as previously described [68]. 
Briefly, a pipette tip was used to pick one isolat- 
ed bacterial colony from a culture agar plate and 
spread it as a thin film on a MTP 384 MALDI- 
TOF target plate [Bruker Daltonics, Leipzig, 
Germany). Twelve distinct deposits from twelve 
isolated colonies were performed for strain 
APIT. Each smear was overlaid with 2 ^iL of ma- 
trix solution [saturated solution of alpha-cyano- 
4-hydroxycinnamic acid) in 50% acetonitrile, 
2.5% tri-fluoracetic acid, and allowed to dry for 
5 minutes. Measurements were performed with 
a Microflex spectrometer [Bruker). Spectra were 
recorded in the positive linear mode for the 
mass range of 2,000 to 20,000 Da [parameter 
settings: ion source 1 [ISI), 20kV; IS2, 18.5 kV; 
lens, 7 kV). A spectrum was obtained after 675 
shots with variable laser power. The time of ac- 
quisition was between 30 seconds and 1 minute 
per spot. The twelve APl^ spectra were import- 



ed into the MALDI BioTyper software [version 
2.0, Bruker) and analyzed by standard pattern 
matching [with default parameter settings) 
against the main spectra of 3,769 bacteria, in- 
cluding 129 spectra from 98 Bacteroides species. 
The method of identification included the m/z 
from 3,000 to 15,000 Da. For every spectrum, a 
maximum of 100 peaks were compared with 
spectra in database. The resulting score enabled 
the identification of tested species, or not: a 
score > 2 with a validly published species ena- 
bled identification at the species level, a score > 
1.7 but < 2 enabled identification at the genus 
level, and a score < 1.7 did not enable any identi- 
fication. No significant MALDI-TOF score was 
obtained for strain APIt against the Bruker da- 
tabase, suggesting that our isolate was not a 
member of a known species. We added the spec- 
trum from strain APl^ to our database [Figure 
4). Finally, the gel view showed the spectral dif- 
ferences with other members of the genus 
Bacteroides [Figure 5). 




Figure 4. Reference mass spectrum from B. timonensis strain APl^. Spectra from 12 individual colonies 
were compared and a reference spectrum was generated. 
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Bacteroides vulgatLS DSM 1447T 
Bacteroides jreolyticus DSM20703T 
Bacteroides uniformis DSM6597T 
Bacteroides timonensis 
Bacteroides thetaiotao micron DSM 2079T 
Bacteroides stercoris DSM 19555T 
Bacteroides salyersiae DSM 18765T 
Bacteroides ovatus DSM 1896T 
Bacteroides nordii DSM 1B764T 
Bacteroides massiliensis DSM 17679T 
Bacteroides intestinalis DSM 17393T 
Bacteroides heparirolyticus 920BB82 
Bacteroides gallinarum DSM 18171T 
Bacteroides frag ilis DSM 9669 
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Bacteroides cellulosilyticus 1BBG9242 
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Figure 5. Gel view comparing B. timonensis strain APl^ to other Bacteroides species. The gel view displays the raw 
spectra of loaded spectrum files as a pseudo-electrophoretic gel. The x-axis records the m/z value. The left y-axis 
displays the running spectrum number originating from subsequent spectra loading. The peak intensity is expressed 
by a grey scale scheme code. The grey scale bar on the right y-axis indicate the relation between the shade of grey 
a peak is displayed with and the peak intensity in arbitrary units. Displayed species are detailed in the left column. 



Genome sequencing information 

Genome project history 

The organism was selected for sequencing on 
the basis of its phylogenetic position and 16S 
rRNA gene sequence similarity to members of 
the genus Bacteroides, and is part of a study of 
the human digestive flora aiming at isolating all 
bacterial species within human feces [1]. It was 



the ninety-ninth genome of a Bacteroides species 
and the first genome of B. timonensis sp. nov. 
The GenBank accession number is 
CBVIOOOOOOOOO and consists of 211 contigs. Ta- 
ble 3 shows the project information and its asso- 
ciation with MIGS version 2.0 comphance [43]. 



Table 3. Project information 



MIGS ID 


Property 


Term 


MIGS-31 


Finishing quality 


High-quality draft 


MIGS-28 


Libraries used 


454 GS paired-end 3-kb library 


MIGS-29 


Sequencing platform 


454 GS FLX Titanium 


MIGS-31. 2 


Fold coverage 


35.76 


MIGS-30 


Assemblers 


gsAssembler 


MIGS-32 


Gene calling method 


PRODIGAL 



Growth conditions and DNA isolation 

B. timonensis sp. nov., strain APl^ [= CSUR P194 
= DSM 26083) was grown on 5% sheep blood- 
enriched Columbia agar (BioMerieux) at 37°C in 
anaerobic atmosphere. Bacteria grown on four 



Petri dishes were harvested and resuspended in 
4xl00nL of TE buffer. Then, 200^1 of this sus- 
pension was diluted in 1ml TE buffer for lysis 
treatment that included a 30- minute incubation 
with 2.5 |ig/|iL lysozyme at 37°C, followed by an 
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overnight incubation with 20 ^g/|iL proteinase 
K at 37°C. Extracted DNA was then purified us- 
ing 3 successive phenol-chloroform extractions 
and ethanol precipitation at -20°C overnight. Af- 
ter centrifugation, the DNA was resuspended in 
160 \iL TE buffer. The yield and concentration 
was measured by the Quant-it Picogreen kit 
(Invitrogen) on the Genios-Tecan fluorometer at 
88.6 ng/nl. 

Genome sequencing and assembly 

Five |ig of DNA was mechanically fragmented on 
Covaris device (KBioScience-LGC Genomics, 
Teddington, UK] using miniTUBE-blue. The DNA 
fragmentation was visualized through an Agilent 
2100 BioAnalyzer on a DNA labchip 7500 with 
an average size of 2.950kb. A 3 kb paired-end li- 
brary was constructed according to the 454 GS 
FLX Titanium paired-end protocol (Roche). Cir- 
cularization and nebulization were performed 
and generated a pattern with a mean size of 513 
bp. After PGR amplification through 17 cycles 
followed by double size selection, the single 
stranded paired-end library was quantified with 
the Quant-it Ribogreen kit [Invitrogen) on the 
Genios Tecan fluorometer at 243 pg/|^L. The li- 
brary concentration equivalence was calculated 
as 8.69 X 108 molecules/ ^iL. The library was 
stored at -20°C until further use. 

The paired-end library was clonally amplified 
with 0.5cpb and Icbp in 8 SV-emPCR reactions 
with the GS Titanium SV emPCR Kit (Lib-L) v2 
[Roche). The yields of the emPCR reactions were 
4.65 and 7.29% respectively, within the recom- 
mended range of 5 to 20% from the Roche pro- 
cedure. Approximately 790,000 beads were 
loaded on a 1/4 region of a GS Titanium 
PicoTiterPlate PTP Kit 70x75 and sequenced 
with the GS Titanium Sequencing Kit XLR70 
[Roche). The run was performed overnight and 
then analyzed on the cluster through the 
gsRunBrowser and Newbler assembler [Roche). 
A total of 802,249 passed filter wells were ob- 
tained and generated 255Mb with a length aver- 
age of 314 bp. These sequences were assembled 
using Newbler [Roche) with 90% identity and 
40bp as overlap. The final assembly identified 
63 scaffolds and 211 large contigs [>l,500bp) 
generating a genome size of 7.13 Mb which cor- 
responds to a coverage of 35.76x genome equiv- 
alent. 

Genome annotation 

Open Reading Frames [ORFs) were predicted us- 
ing Prodigal [69] with default parameters. How- 
ever, the predicted ORFs were excluded if they 
spanned a sequencing gap region. The predicted 
bacterial protein sequences were searched 



against the GenBank [70] and Clusters of 
Orthologous Groups [COG) databases using 
BLASTP. The tRNAs and rRNAs were predicted 
using the tRNAScan-SE [71] and RNAmmer [72] 
tools, respectively. Signal peptides and numbers 
of transmembrane helices were predicted using 
SignalP [73] and TMHMM [74], respectively. 
Mobile genetic elements were predicted using 
PHAST [75] and RAST [76]. ORFans were identi- 
fied if their BLASTP E- value was lower than le- 
03 for alignment length greater than 80 amino 
acids. If alignment lengths were smaller than 80 
amino acids, we used an E-value of le-05. Such 
parameter thresholds have already been used in 
previous works to define ORFans. Artemis [77] 
and DNA Plotter [78] were used for data man- 
agement and visualization of genomic features, 
respectively. The Mauve alignment tool [version 
2.3.1) was used for multiple genomic sequence 
alignment [79]. 

To estimate the mean level of nucleotide se- 
quence similarity at the genome level between B. 
timonensis and 9 other members of the genus 
Bacteroides [Table 6), we used the Average Ge- 
nomic Identity Of gene Sequences [AGIOS) in- 
house software [8]. Briefly, this software uses 
the Proteinortho software [80] for the pairwise 
detection of orthologous proteins between ge- 
nomes, then retrieves the corresponding genes 
and determines the mean percentage of nucleo- 
tide sequence identity among orthologous ORFs 
using the Needleman-Wunsch global alignment 
algorithm. B. timonensis strain APIt was com- 
pared to B. intestinalis strain DSM 17393 
[GenBank accession number 

NZ_ABJLOOOOOOOO), B. ceUnIosilyticus strain 
DSM 14838 [NZ.ACCHOOOOOOOO), B fragilis 
strain YCH46 [NC_006347), B. vulgatus strain 
ATCC 8482 [NC_009614), B. thetaiotaomicron 
strain VPI-5482 [NC_004663), B. salanitronis 
strain DSM 18170 [NC_015164), B. helcogenes 
strain P36-108 [NC_014933), B. flnegoldii strain 
DSM 17565 [NZ.ABXIOOOOOOOO) and B. 
uniformis strain ATCC 8492 [AAYHOOOOOOOO). 

Genome properties 

The genome is 7,130,768 bp long [1 chromo- 
some, but no plasmid) with a 43.3% G+C content 
[Figure 6 and Table 4). Of the 5,845 predicted 
genes, 5,786 were protein-coding genes and 59 
were RNAs, including 1 complete rRNA operon. 
A total of 3,111 genes [53.22%) were assigned a 
putative function and 3,283 genes were identi- 
fied as ORFans [56.16%). Strain APIt possesses 
a variety of mobile genetic elements. These in- 
clude 6 prophages of 13.70, 14.60, 10.51, 8.18, 
9.91 and 12.79 Kb, respectively) and 91 trans- 
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posable elements belonging to 18 transposon 
families that include the putative mobilization 
protein BF0133, the putative conjugative trans- 
poson mobihzation protein BF0132, the hypo- 
thetical protein clustered with conjugative 
transposons BF0131, TraA-CTn, TraB-CTn,TraD- 
CTn,TraE-CTn,TraF-CTn,TraG-CTn,TraH- 

7100000 



I I 




3550000 



Figure 6. Graphical circular map of the chromosome. From the outside in: open reading frames oriented 
in the forward (colored by COG categories) direction, open reading frames oriented in the reverse (col- 
ored by COG categories) direction, RNA operon (red), and tRNAs (green), GC content plot, and GC skew 
(purple: negative values, olive: positive values). 



Table 4. Nucleotide content and gene count 


levels of the genome 




Attribute 


Value 


% of totaP 


Genome size (bp) 


7,130,768 




DNA coding region (bp) 


6,434,142 


90.23 


DNA G+C content (bp) 


3,087,622 


43.30 


Number of replicons 


1 




Extra chromosomal element 


0 




Total genes 


5,845 


100 


RNA genes 


59 


1.01 


Protein-coding genes 


5,786 


98.99 


Genes with function prediction 


3,111 


53.22 


Genes assigned to COGs 


2,820 


48.24 


Genes with peptide signals 


435 


7.44 


Genes with transmembrane helices 


456 


7.80 


" The total is based on either the size of the 


genome in base pairs or the total 


number of pro- 


tein-coding genes in the annotated genome 







CTn,TraI-CTn,TraJ-CTn,TraK-CTn,TraL- 
CTn,TraM-CTn,TraN-CTn,TraO-CTn and TraQ- 
CTn. The properties and statistics of the genome 
are summarized in Tables 4 and 5. The distribu- 
tion of genes into COGs functional categories is 
presented in Table 5. 
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Table 5. Number of genes associated with the 25 


general COG functional categories 


Code 


Value 


%age" 


Description 


J 


156 


2.66 


Translation 


A 


0 


0 


RNA processing and modification 


K 


234 


4.00 


Transcription 


L 


200 


3.42 


Replication, recombination and repair 


B 


0 


0 


Chromatin structure and dynamics 


D 


27 


0.46 


Cell cycle control, mitosis and meiosis 


Y 


0 


0 


Nuclear structure 


V 


107 


1.83 


Defense mechanisms 


T 


240 


4.22 


Signal transduction mechanisms 


M 


361 


6.17 


Cell wall/membrane biogenesis 


N 


5 


0.08 


Cell motility 


z 


0 


0 


Cytoskeleton 


w 


0 


0 


nxLi dcei 1 u idF siruLLures 


u 


65 


1.11 


II 1 Lfdcei lu lar Lraiiiciviiig aiiu becreuuii 








Posttranslational modification, protein turnover, chaper- 


o 


89 


1.52 


ones 


c 


168 


2.87 


Energy production and conversion 


G 


369 


6.31 


Carbohydrate transport and metabolism 


E 


212 


3.62 


Amino acid transport and metabolism 


F 


73 


1.25 


Nucleotide transport and metabolism 


H 


130 


2.22 


Coenzyme transport and metabolism 


1 


87 


1.48 


Lipiu transport anu meiaooiism 


P 


202 


3.42 


IIIUIgdlllL lOII LIdllSUUI I dllU IlieidUUl ISIII 








Secondary metabolites biosynthesis, transport and ca- 


Q 


47 


0.80 


tabolism 


R 


518 


8.86 


General function prediction only 


S 


197 


3.37 


Function unknown 




2966 


51 .26 


Not in COGs 



^ The total is based on the total number of protein-coding genes in the annotated genome. 



Genome comparison with other 
Bacteroides genomes 

Here, we compare the genome of B. timonensis 
with those of B. intestinalis, DSM 17393, B. 
cellulosilyticus DSM 14838, B. fragilis YCH46, B. 
vulgatus ATCC 8482, B. thetaiotaomicron VPI- 
5482, B. salanitronis DSM 18170, B. helcogenes P 
36-108, B.finegoldii DSM 17565 and B. uniformis 
ATCC 8492. The draft genome of B. timonensis 
(7.13Mb) is larger than all other studied ge- 
nomes (Table 6A]. It also exhibits a higher G+C 
content than all other genomes exceptB. 
salanitronis, B. helcogenes and B. uniformis (43.3, 
46.4, 44.7 and 46.4%, respectively). B. 
timonensis has a higher gene content (5,786) 
than any other compared genome. The distribu- 
tion of genes into COG categories was similar in 
all 10 compared genomes except in the N cate- 



gory (cell motility) for which B. fragilis, B. 
vulgatus, B. salanitronis, B. helcogenes and B. 
uniformis were underrepresented (Figure 7). In 
addition, B. timonensis shared 2,956, 3,081, 
2,159, 2,099, 2,379, 1,721, 2,001, 2,039 and 
2,268 orthologous genes with B. intestinalis, B. 
cellulosilyticus, B. fragilis, B. vulgatus, B. 
thetaiotaomicron, B. salanitronis, B. helcogenes, 
B. finegoldii and B. uniformis, respectively. 
Among compared genomes except B. timonensis, 
AGIOS values ranged from 70.16 between B. 
salitronis and B. cellulosilyticus to 88.16% be- 
tween B. intestinalis and B. cellulosilyticus. When 
B. timonensis was compared to other species, 
AGIOS values ranged from 70.29 with B. 
salitronis to 93.61% with B. cellulosilyticus (Ta- 
ble 6B). 
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Table 6A. Genomic 


comparison of B. timonensis with 9 other Bacteroides species'*'. 






Species 






Strain 




Genome acces- 
sion 
number 


Genome size 
(Mb) 


G+C content 


B. timonensis 
B. intestinalis 
B. cellulosilyticus 




API 
DSM 17393 
DSM 14838 


CBVI01 0000000 
NZ_ABJLOOOOOO 

00 

NZ_ACCHOOOOO 
000 




7.13 
6.05 
6.87 




43.3 
42.8 
42.7 


B. fragilis 






YCH46 




NC 006347 




5.28 




43.2 


B. vulgatus 






ATCC 8482 


NC 009614 




5.16 




42.2 


B. thetaiotaomicron 




VPI-5482 


NC_004663 




6.26 




42.8 


B. salanitronis 




DSM 18170 


NC_015164 




4.24 




46.4 


B. helcogenes 
B. finegoldii 




P 36-108 
DSM 17565 


NC_014933 

NZ_ABX1000000 

00 




4.0 
4.89 




44.7 
42.9 


B. uniformis 




ATCC 8492 


AAYHOOOOOOOO 




4.72 




46.4 


^Species, Strain, GenBank accession number, genome size and G+C content of 


all compared genomes. 


Table 6B. Genomic 


; compar 


ison of 




9 oth 


er 


species'. 








B. tint 


B. int 


fi. eel 


fi. fra 


fi. vul 


fi. the 


fi. sal 


fi. hel 


fi. fin 


fi. uni 


B. tim 


5,786 


2,956 


3,081 


2,159 


2,099 


2,379 


1,721 


2,001 


2,039 


2,268 


B. int 


87.73 


4,911 


2,967 


2,085 


2,036 


2,361 


1,667 


1,963 


2,066 


2,278 


B. cell 


93.61 


88.16 


5,719 


2,130 


2,078 


2,380 


1,655 


1,990 


2,017 


2,231 


B. fra 


73.76 


74.43 


73.92 


4,184 


1,927 


2,174 


1,517 


1,893 


1,880 


1,995 


B. vul 


7^m 


71.74 


71.48 


71.87 


4,066 


2,100 


1,638 


1,743 


1,859 


1,898 


B. the 


73.99 


74.65 


73.87 


75.42 


72.21 


4,778 


1,601 


1,891 


2,191 


2,039 


B. sal 


70.29 


70.65 


70.16 


70.35 


72.18 


70.50 


3,553 


1,466 


1,580 


1,584 


B. hel 


76.40 


76.51 


76.41 


74.15 


71.62 


73.64 


70.68 


3,244 


1,703 


1,930 


B. fin 


74.28 


75.01 


74.45 


75.72 


72.22 


81.24 


70.77 


73.99 


4,485 


1,920 


B. uni 


77.08 


76.83 


76.80 


74.25 


72.45 


74.36 


71.32 


79.37 


74.77 


4,663 



t numbers of orthologous proteins shared between genomes (above diagonal), AGIOS values (below di- 
agonal) and numbers of proteins per genome (bold numbers). 



B. tim = B. timonensis, B. int = B. intestinalis, B. eel = B. cellulosilyticus, B. fra = B. fragilis, B. vul = B. 
vulgatus, B. the = B. thetaiotaomicron, B. sal = B. salanitronis, B. hel = B. helcogenes, B. uni = B. 
uniformis, B. fin = B. finegoldii. 
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I B.uni 
I B.fin 

B.hel 
I B.sal 

B.the 
I B.vul 
I B.fra 

B.cel 
I B.int 
I B.tim 



JAKLBDYVTMNZWUOCGEFH I PQRS 

COG categories 

Figure 7. Distribution of predicted genes of B. timonensis and 9 other Bacteroides species into COG cat- 
egories. B. uni = B. uniformis, B. fin = B. finegoldii, B. hel = B. helcogenes, B. sal = B. salanitronis, B. the 
= B. thetaiotaomicron, B. vul = B. vulgatus, B. fra = B. fragilis, B. eel = B. cellulosilyticus, B. int = B. 
intestinalis, B. tim = B. timonensis. 



Conclusion 

On the basis of phenotypic, phylogenetic and ge- 
nomic analyses [taxono-genomics), we formally 
propose the creation of Bacteroides timonensis 
sp. nov. that contains strain APIt. This strain 
was isolated from the fecal flora of a 21-year-old 
woman who suffered from severe anorexia ner- 
vosa. 

Description of B. timonensis sp. nov. 

Bacteroides timonensis [tim.o.nen'sis. L. masc. 
adj. timonensis, of Timone, the name of the hos- 
pital where strain APl^ was first cultivated). 

Colonies are translucent and 0.3 mm in diameter 
on blood-enriched Columbia agar. Cells are rod- 
shaped with a mean diameter of 0.88 ^im. Opti- 
mal growth is achieved anaerobically, although 
the strain is able to grow under microaerophilic 
conditions, and weakly with 5% CO2. Growth oc- 
curs between 25°C and 37°C, with optimal 
growth at 37°C. Cells stain Gram-negative and 
are not motile. Positive reactions for catalase, 
arginine dihydrolase, a-galactosidase, (B- 
galactosidase, a-glucosidase, (B-glucosidase, a- 
arabinosidase, N-acetyl-(3-glucosaminidase, glu- 
tamic acid decarboxylase, a-fucosidase, nitrate 
reduction, indole production, alkaline phospha- 
tase, proline arylamidase, leucyl glycine 



arylamidase, alanine arylamidase, glutamyl glu- 
tamic acid arylamidase, and fermentation of 
mannose and raffinose. 

Weak activities are observed for glycine 
arylamidase and serine arylamidase. Negative 
reactions are obtained for urease, (B- 
galctosidase-6-phosphatase, (B-glucuronidase, 
arginine arylamidase, phenylalanine 
arylamidase, leucine arylamidase, pyroglutamic 
acid arylamidase, tyrosine arylamidase and 
histidine arylamidase. Using an API 50CH strip 
[Biomerieux), strain APIt is asaccharo lytic. Cells 
are susceptible to susceptible to amoxicillin- 
clavulanate, ceftriaxone, imipenem, trime- 
thoprim-sulfamethoxazole, metronidazole and 
doxycycline but resistant to amoxicillin, 
vancomycin and gentamicin. 

The 16S rRNA and genome sequences are depos- 
ited in GenBank under accession numbers 
JX041639 and CBVIOOOOOOOOO, respectively. 
The G+C content of the genome is 43.3%. The 
habitat of the organism is the digestive tract. The 
type strain APIt (= CSUR P194 = DSMZ 26083) 
was isolated from the fecal flora of a French 
Caucasoid female who suffered from a severe re- 
strictive form of anorexia nervosa. This strain 
has been found in Marseille, France. 
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